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ABSTRACT 

We use a theoretical model to predict the clustering properties of galaxy clus- 
ters. Our technique accounts for past light-cone effects on the observed cluster- 
ing and follows the non-linear evolution in redshift of the underlying dark mat- 
ter correlation function and cluster bias factor. A linear treatment of redshift- 
space distortions is also included. We perform a maximum-likelihood analysis 
by comparing the theoretical predictions to a set of observational data, both 
in the optical (two different subsamples of the APM catalogue and the EDCC 
catalogue) and X-ray band (RASSl Bright Sample, BCS, XBACs, REFLEX). 
In the framework of cold dark matter models, we compute the constraints on 
cosmological parameters, such as the matter density f^om, the cosmological 
constant f2oA, the power-spectrum shape parameter T and normalisation erg. 
Our results show that X-ray data are more powerful than optical ones, allow- 
ing smaller regions in the parameter space. If we fix T and cts to the values 
suggested by different observational datasets, we obtain strong constraints on 
the matter density parameter: fJom < 0.5 and 0.2 < fJom < 0.35, for the op- 
tical and X-ray data, respectively. Allowing the shape parameter to vary, we 
find that the clustering properties of clusters are almost independent of the 
matter density parameter and of the presence of a cosmological constant, while 
they appear to be strongly dependent on the shape parameter. Using the X- 
ray data only, we obtain F ~ 0.1 and 0.4^ crg^ 1.1 for the Einstein-de Sitter 
model, while 0.14^ F^ 0.22 and 0.6^ (Tg^ 1.3 for open and flat models with 
^Om = 0.3. Finally, we use our model to make predictions on the correlation 
length of galaxy clusters expected in future surveys. In particular, we show the 
results for an optical catalogue with characteristics similar to the EIS project 
and for a very deep X-ray catalogue with the characteristics of the XMM/LSS 
survey. We find that clusters at high redshifts are expected to have larger a 
correlation length than local ones. 

Key words: cosmology: theory - galaxies: clusters - large-scale structure of 
Universe - X-rays: galaxies - dark matter 



1 INTRODUCTION 

Clusters of galaxies are the largest collapsed objects 
with masses dominated by dark matter. Their proper- 
ties are largely determined by the gravitational collapse 
in the cosmological density field, and so their distribu- 
tion in the Universe depends both on cosmology and on 
the initial density perturbations. Therefore, the study 
of the clustering of clusters may provide important di- 



agnostics for models of structure formation. One impor- 
tant advantage of using clusters to study the large-scale 
structure is that their formation in the cosmic density 
field is relatively easy to understand. Indeed, modern 
cosmological simulations show that clusters of galaxies 
may be identified as the most massive dark haloes pro- 
duced by gravitational collapse, and so their clustering 
properties are relatively easy to interpret theoretically. 
In fact, based on an extension of the Press-Schechter 
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formalism, one can obtain an analytical model in which 
the dependence of the two-point correlation function of 
dark haloes on cosmology and initial power spectrum 
can be explicitly seen (Mo & White 1996; Catelan et al. 
1998; Sheth, Mo & Tormen 2001). 

It has been known for some time that the strong 
observed correlation of clusters on large scales is very 
difficult to reconcile with the standard cold dark mat- 
ter (CDM) model (White et al. 1987; Dahon et al. 1992; 
Jing et al. 1993; Mo, Peacock & Xia 1993; Dalton et al. 
1994; Borgani, Coles & Moscardini 1994; Croft & Efs- 
tathiou 1994; Borgani et al. 1995). More recently, anal- 
yses on the constraints by the cluster-cluster correlation 
function on current theoretical models of structure for- 
mation have been carried out by a number of authors 
(Eke et al. 1996; Mo, Jing & White 1996; Moscardini et 
al. 2000a,b; Robinson 2000; Colberg et al. 2000). 

In this paper, we use current available observational 
data on the cluster two-point correlation function to 
constrain theoretical models of structure formation. We 
use analytical tools well tested by numerical simulations 
to make predictions for a large grid of models. Our anal- 
ysis follows closely that of Mo et al. (1996) and Robin- 
son (2000). However, there are several important differ- 
ences. First, we consider new datasets; in particular, we 
use data from recent X-ray surveys (see also Moscardini 
et al. 2000b) . Second, we use improved theoretical mod- 
els (Matarrese et al. 1997). In particular we pay much 
attention to the redshift evolution in the clustering of 
clusters and to light-cone and selection effects. Third, 
we present theoretical predictions for some future sur- 
veys. 

The plan of the paper is as follows. In Section 2 
we list the observational cluster samples used in the 
following analysis. Section 3 is devoted to the presen- 
tation of our theoretical model to estimate the correla- 
tion of galaxy clusters both in the optical and in the 
X-ray bands. In Section 4 we present the results of 
a maximum-likelihood analysis performed on existing 
data and we discuss the constraints on the cosmological 
parameters. In Section 5 we show the predictions of the 
correlation length for possible future surveys. Conclu- 
sions are drawn in Section 6. 



• To construct the APM cluster survey, an auto- 
mated procedure has been used to select from the an- 
gular APM galaxy survey the clusters whose redshifts 
were later measured. An early sample, containing 364 
clusters with richness TZ > 50, called Sample B (Dal- 
ton et al. 1994), has been created by fixing the limiting 
magnitude of the galaxy catalogue to bj = 20.5. A new, 
richer sample was then obtained by Croft et al. (1997) 
by extending the limit to 6j = 21.0 and considering only 
clusters with richness 7?, > 80. This catalogue (called 
Sample C) contains 165 objects with redshift cz < 55000 
km s~^. Even if the clustering study made by Croft et 
al. (1997), using a maximum-likelihood technique, con- 
sidered six different subsamples with decreasing cluster 
density, in our statistical analysis we will prefer to use 
the results only for the two whole samples B and C. 
This is done to reduce the possible problems coming 
from the assumption of independence of the different 
samples. This assumption is required by the way we 
write the likelihood (see equation H). Our final results 
will not be affected by this choice, as shown by the fact 
that the confidence levels on the cosmological parame- 
ters obtained using sample B only are almost indistin- 
guishable from those coming from the complete analysis 
of the optical datasets. The corresponding mean intra- 
cluster separations for samples B and C are 30 and 57 
h~^ Mpc, respectively {h is the value of the local Hubble 
constant Ho in units of 100 km s^^ Mpc^^); the values 
of the correlation lengths ro, of the slope 7 of the corre- 
lation function, and their la errorbars, obtained from a 
maximum likelihood analysis, are reported in Table m. 

• The EDCC is a machine-based, objectively selected 
sample of galaxy clusters consisting of 737 objects of all 
richness, over 0.5 sr of sky centred on the South Galactic 
Pole (Lumsden et al. 1992). A subsample of 79 clusters, 
with at least 22 galaxies inside a radius corresponding to 
1 h~^ Mpc, with magnitudes between the limits 7713 and 
ms + 2 (see Abell 1958 for the definitions) has been spec- 
troscopically confirmed. The resulting catalogue, with a 
corresponding mean intracluster separation dc — AGh"^ 
Mpc, has been used by Nichol et al. (1992) to estimate 
the clustering properties. The values of ro and 7 ob- 
tained from a least-squares fit are reported in Table hi 



2 CLUSTER SAMPLES 

In this section we present the observational data on the 
galaxy cluster correlation length that we will use in the 
following analysis to constrain the cosmological param- 
eters. A summary of these data is presented in Table 





2.1 Optical data 

As far as the optical band is concerned, we will use for 
our study the observational results from the analysis by 
Croft et al. (1997) on the APM cluster redshift survey 
and by Nichol et al. (1992) on the Edinburgh- Durham 
cluster catalogue (EDCC). 



2.2 X-ray data 

In the following analysis we will consider four different 
catalogues of galaxy clusters selected in the X-ray band: 
the ROSAT All-Sky Survey 1 (RASSl) Bright Sample, 
the ROSAT Brightest Cluster Sample (BCS), the X- 
ray brightest Abell-type cluster sample (XBACs), the 
ROSAT -E SO Flux-limited X-ray sample (REFLEX). 
Here we list some of their characteristics; more details 
can be found in the original papers. 

• The RASSl Bright Sample (De Grandi et al. 1999a) 
contains 130 clusters of galaxies selected from the 
ROSAT All-Sky Survey (RASS) data. The catalogue 
has an effective flux limit in the (0.5 - 2.0 keV) band 
between 3.05 and 4 x 10~^^ erg cm~^ s~^ over the se- 
lected area which covers a region of approximately 2.5 
sr within the Southern Galactic Cap, i.e. S < 2.5° and 
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Table 1. The clustering data (in the optical and X-ray bands) used in the likelihood analysis. Column 1: catalogue name. 
Column 2: characteristics of the catalogue: for the optical catalogues, the richness TZ and the mean intracluster distance dc; for 
the X-ray catalogues, the limiting luminosity Lx or the limiting flux S'limi unless the whole catalogue is analyzed. Column 3: 
reference for the clustering analysis. Column 4: number of clusters in the catalogue ric- Columns 5, 6 and 7: correlation length 
tq (in h^^ Mpc), the slope of the correlation function ■y and the corresponding confidence levels of the quoted crrorbars. 



Catalogue 


Characteristics 


Reference 


ric 


ro (/i-i Mpc) 


7 


Errorbars 


Optical band: 














APM Sample B 


TZ > 50, dc = 30 h-'^ Mpc 


Croft et al. (1997) 


364 


,. 2+0-4 
-^^•^-0.6 


r, -|o + 0.09 

^•-^■^-0.06 


la- 


APM Sample C 


7^ > 80, dc = 57 h-^ Mpc 


idem 


110 


IS 4+2-2 


1 y+0.3 
-^■'-0.3 


id 


EDCC 


dc = 46 h-^ Mpc 


Nichol et al. (1992) 


79 


16.4 ±4.0 


2.1 ±0.3 


la 


X-ray band: 














RASSl Bright Sample 


whole catalogue 


Moscardini et al. (2000a) 


130 


21.5ti1 


2 11 +0-53 
^■11-0.56 


2a 


BCS 


Lx > 0.24 X W^^h-'^ erg s^^ 


Lee & Park (1999) 


33 


33.0«-J 


-I 09+0.49 


la 


XBACs 


whole catalogue 


Borgani et al. (1999) 


203 


26.0 ±4.5 


1 98+°-^'' 


2a 


REFLEX 


Sn^ = 3 X 10-12 erg cm-2 s-i 


Collins et al. (2000) 


449 


18.8 ±0.9 


1 SQ+0.15 

-^•"•^-o.os 


la 



b < —20°. For our theoretical predictions we will use the 
exact sky map of the sample which is presented in Figure 
2 of De Grandi et al. (1999a). The redshift distribution 
has a tail up to z ~ 0.3 but the majority of the clus- 
ters have z < 0.1. Moscardini et al. (2000a) found that 
the two-point correlation function of the whole sample 
is well fitted by a power-law with ro — 21.5^\''^h~^ Mpc 



and 7 = 2.11 



+0.53 



(95.4 per cent confidence level with 



one fitting parameter). 

• The BCS catalogue (Ebeling et al. 1998) is an X- 
ray selected, flux-limited sample of 201 galaxy clusters 
with z < 0.3 drawn from the RASS data in the north- 
ern hemisphere {S > 0°) and at high Galactic latitude 
(I&//I > 20°). The limiting flux is Sun, = 4.45 x 10"^^ 
erg cm~^ s~i in the (0.1-2.4) keV band. Since the sky- 
coverage flskyiS) of BCS is not available, we will use 
^sky(>S') = const ~ 4.13 steradians for fluxes larger 
than S'lim- Lee & Park (1999) analyzed the cluster- 
ing properties of this catalogue using four different 
volume-limited subsamples. For simplicity, in the fol- 
lowing analysis we will consider only the catalogue with 
Lx > 0.24 X lO'^'^h^'^ erg s"^ and a limiting redshift 
oi z = 0.07. The correlation function in this case is fit- 
ted by a power-law with ro — SS.Olgg/i"^ Mpc and 
7 — 1.82t.Q''^Q (68.3 per cent confidence level with one 
fitting parameter). We checked that the inclusion of the 
deeper subsamples produced very small changes in our 
likelihood analysis. 

• The XBACs catalogue (Ebeling et al. 1996) is an 
all-sky X-ray sample of 242 Abell galaxy clusters ex- 
tracted from the RASS data. Being optically selected, 
it is not a complete flux-limited catalogue. The sam- 
ple covers high Galactic latitudes (|bi/| > 20°). The 
adopted limiting flux is Sum = 5 x 10^^^ erg cm^^ 
s^i in the (0.1-2.4) keV band. Also in this case, since 
the actual sky coverage is not available, we will adopt 
fisky('S') — const ~ 8.27 steradians for fluxes larger than 
Slim. The aforementioned selection effects produce a lu- 
minosity function for XBACs which is much lower in the 
faint part than that obtained from other catalogues (e.g. 
Ebeling et al. 1997; Rosati et al. 1998; De Grandi et al. 



1999b). We take into account this incompleteness in our 
model following the same method described in Moscar- 
dini et al. (2000b). The clustering properties of this cat- 
alogue have been studied by different authors. Abadi, 
Lambas & Muriel (1998) found for the whole catalogue 
a correlation length ro = 21.1123^^^ Mpc (la errorbar) 
and a slope 7 — —1.92. Borgani, Plionis & Kolokotro- 
nis (1999) analyzed the same sample finding a somewhat 
larger correlation amplitude: ro = 26.0'^^}^h~^ Mpc and 
7 = 1.981^0 53 (95.4 per cent confidence level with one 
fitting parameter). The difference between these two es- 
timates is probably due to a different assumption on 
the errors. Since the maximum-likelihood approach is 
more robust than the quasi-Poisson assumption made 
by Abadi et al. (1998), we prefer to use in the following 
analysis the results obtained by Borgani et al. (1999). 

• The REFLEX survey (Bohringer et al. 1998) is 
a large sample of optically confirmed X-ray clusters 
selected from RASS. The sample includes 452 X-ray 
selected clusters in the southern hemisphere, at high 
Galactic latitude (|b//| > 20°). For our computations, 
we use the actual sky coverage given in Figure 1 of 
Collins et al. (2000). Using a catalogue with a limiting 
fiux of 3 X 10"^^ erg cm~^ s^^ (defined in the 0.1 - 2.4 
keV energy band), where the sky coverage falls to 97.4 
per cent of the whole surveyed region (4.24 steradians), 
Collins et al. (2000) found that the two-point correlation 
function is fitted by a power-law with ro — 18.8±0.9/i~^ 
Mpc and 7 = 1.83lo 08 (68.3 per cent confidence level 
from a maximum- likelihood analysis). As for the APM 
subsamples at higher richness, in our analysis we prefer 
not to use the clustering results for the REFLEX sub- 
samples at higher luminosities, also discussed in Collins 
et al. (2000). 



3 THEORETICAL MODELS FOR THE 
CORRELATION FUNCTION 

Our theoretical predictions for the spatial two-point cor- 
relation function of galaxy clusters in different cosmo- 
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logical models have been here obtained by using an up- 
dated version of the method presented in Moscardini 
et al. (2000a,b), where the application was limited to 
X-ray selected clusters. Here we will give only a short 
description of the method and we refer to those papers 
for a more detailed discussion. 



3.1 Clustering in the past-light cone 

Matarrese et al. (1997; see also Moscardini et al. 1998 
and Yamamoto & Suto 1999 and references therein) de- 
veloped an algorithm to describe the clustering in our 
past light-cone taking into account both the non-linear 
dynamics of the dark matter distribution and the red- 
shift evolution of the bias factor. The final expression 
for the observed spatial correlation function ^obs in a 
given redshift interval Z is 



Cobs r) = ^ ; — -2 

[S,dz,M{z,)Y 



(1) 



where N{z) = M{z)/r{z) and N{z) is the actual red- 
shift distribution of the catalogue. In the previous for- 
mula £,ohi{r,zi,Z2) represents the correlation function 
of pairs of objects at redshifts z\ and 22 with comoving 
separation r. An accurate approximation for it is given 
by 



£,ohi{r,Zl,Z2) ~ &cff(2l)fecff(z2)Cm(r, Zavo) 



(2) 



where ^m is the dark matter covariance function and 
Zavc is a suitably defined intermediate redshift. 

Another important ingredient entering the previous 
equation is the effective bias b^s which can be expressed 
as a weighted average of the 'monochromatic' bias factor 
b{M, z) of objects of some given intrinsic property M 
(like mass, luminosity, etc): 

b,«(z)=U{zy^ I dhiAf b{M',z)Af{z,M'), (3) 

J M 

where J^{z, M) is the number of objects actually present 
in the catalogue with redshift in the range z, z + dz and 
M in the range M, M + dM, whose integral over In M 
isN{z). 

As galaxy clusters are expected to form by the hi- 
erarchical merging of smaller mass units, one can fully 
characterize their properties by the mass M of their 
hosting dark matter haloes at each redshift z. One can 
then estimate their comoving mass function n{z, M) by 
the Press-Schechter (1974) formula and adopt the Mo & 
White (1996) relation for the monochromatic bias. Ac- 
tually, we use the relations recently introduced by Sheth 
& Tormen (1999) and Sheth, Mo & Tormen (2001), 
which have been shown to produce a more accurate fit 
of the distribution of the halo populations in numerical 
simulations (Jenkins et al. 2001). 

To complete our model we need a technique to com- 
pute the redshift evolution of the dark matter covariance 
function fm. For this, we use the fitting formula by Pea- 
cock & Dodds (1996), which allows to evolve ^m into 
the fully non-linear regime. 



3.2 From the catalogue characteristics to the 
halo mass 

In order to predict the abundance and clustering of 
galaxy clusters in the different catalogues described 
above (both in the optical and X-ray bands) we need 
to relate the sample characteristics to a corresponding 
halo mass at each redshift. 

For the optical catalogues, the various samples of 
APM and EDCC considered by Croft et al. (1997) and 
Nichol et al. (1992) respectively, are characterized by a 
different cluster number density (see Table hi). In this 
case, for each cluster subsample, we fix the minimum 
mass (which in our model represents the lower limit of 
the integral in equation H) of the hosting dark matter 
haloes in such a way that the comoving cumulative mass 
function reproduces the number density of clusters in 
the range of redshift sampled by the original catalogues. 

For X-ray selected clusters, we use a different ap- 
proach (see Moscardini et al. 2000a,b and Suto et al. 
2000). In this case it is quite easy to relate the limit- 
ing flux (and/or eventually the limiting luminosity, as 
required for instance by the BCS subsample that we 
consider) to the minimum halo mass. In fact the flux 

5 in a given band corresponds to an X-ray luminosity 
Lx ~ indj^S in the same band, where dL is the luminos- 
ity distance. The quantity Lx can be converted into the 
total luminosity Lboi by performing band and bolomet- 
ric corrections (we assume an overall ICM metallicity 
of 0.3 times solar). Local observations suggest that the 
cluster bolometric luminosity is related to the tempera- 
ture by a simple relation: T — A iboii where the temper- 
ature is expressed in keV and Lboi is in units of 10*'* /i~^ 
erg s~^. We assume A — 4.2 and B — 1/3, which are a 
good representation of the data with T^ 1 keV (e.g. 
Markevitch 1998 and references therein). Below this 
temperature (i.e. for galaxy groups) the Lboi — T re- 
lation has much steeper a slope: for this reason we set a 
minimum temperature at T = 1 keV. Recent analyses of 
cluster temperature data at higher redshifts (Mushotzky 

6 Scharf 1997; Donahue et al. 1999) are consistent with 
no evolution in the Lboi — T relation out to 2 ~ 0.4. 
Therefore, we can safely assume that the previous re- 
lation holds true in the redshift range sampled by the 
considered catalogues. Moscardini et al. (2000b), who 
allowed a moderate redshift evolution of the Lboi — T 
to reproduce the observed log A''-log S relation, showed 
that the clustering properties are only slightly sensitive 
to this assumption. To convert the cluster temperature 
into the mass of the hosting dark matter halo, we assume 
virial isothermal gas distribution and spherical collapse: 
T oc M^^^E^^^{z)A^ir{zy^^, where Avir represents the 
mean density of the virialized halo in units of the critical 
density at that redshift and E{z) is the ratio between 
the value of the Hubble constant at redshift z and to- 
day. Once the relation between observed flux and halo 
mass at each redshift is established we have to account 
for the catalogue sky coverage f2sky (S) (when available) 
to predict the redshift distribution. 
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3.3 Redshift distortions 

We include the effect of redshift-space distortions using 
linear theory and the distant-observer approximation 
(Kaiser 1987). Under these assumptions the enhance- 
ment of the redshift-space averaged power spectrum is 
approximately 1-1-2/3/3-1-/3^/5, where /3 ~ Qo^/b^ff. We 
find that the total effect of redshift distortions on the 
value of the correlation length is always smaller than 10 
per cent. 



4 RESULTS 

4.1 Cosmological models 

In what follows we use a maximum-likelihood analysis 
to constrain the main parameters defining a cosmolog- 
ical model. We will consider a set of structure forma- 
tion models all belonging to the general class of cold 
dark matter (CDM) ones, for which the linear power- 
spectrum reads Piin(fc,0) oc k"T^{k), with T(fc) the 
CDM transfer function (Bardeen et al. 1986) 

T(.) ^ '-^\V:'''^ X [1-^3.89.-, 



2.34q 



(4) 



where q = k/hV. The shape parameter F depends on the 
Hubble parameter h, on the matter density fiom and on 
the baryon density f2ob (Sugiyama 1995): 



r = r2om/iexp(-f2ob - vVO^ ^ob/f2om) 



(5) 



We fix the spectral index n to unity and we allow 
r to vary in the range 0.05-0.5, while nom ranges from 
0.1 to 1 in the framework of both open and flat models, 
with a cosmological constant contributing to the total 
density with f7oA = 1 — f^Om. Finally, we use different 
normalisations of the primordial power-spectrum, pa- 
rameterized by (78 (the r.m.s. fluctuation amplitude in 
a sphere of 8/i~^ Mpc) in the range 0.2 < erg < 2. 

In summary, the cosmological models we consider 
are defined by four parameters: f2om, f2oA, F and erg. 



4.2 Maximum likelihood analysis 

Confidence levels for the cosmological parameters are 
obtained through a maximum likelihood analysis. The 
likelihood is £ oc exp(— x^/2), where 



v^'" [ro{i) - ro(i;riom,f2oA,F,CTg)]^ 
1=1 



Jro 



(i) 



(6) 



The sum runs over the observational dataset described 
in Section 2, i.e. A'data ~ 3 and A'^data = 4 for the optical 
and X-ray bands, respectively. The quantities ro{i) and 
arg (i) represent the values of the correlation length and 
its la errorbar for each catalogue, as reported in Table 
m ro{i; f^om, f^oA, T, erg) is the corresponding theoretical 
prediction obtained with a given choice of cosmological 
parameters. Note that the observed correlation lengths 
reported in Table hi have been obtained using different 
fitting methods (least-squares or maximum-likelihood 



analyses) and over different spatial scales in ways which 
could be inconsistent with each other. Moreover, the 
quoted values of ro come from fits of the data where also 
the values of the slope are marginalized and it is known 
that there is significant covariance between ro and 7. 
In spite of this, we prefer not to estimate ro through 
the standard power-law fit procedure with two param- 
eters, because in the case of the maximum-likelihood 
analyses the comparison would not be appropriate. For 
these reasons we deflne the predicted correlation length 
as the distance where the two-point correlation function 
is unity. However, in all cases we also checked that the 
values of the slope 7 obtained from the power-law flt 
are well inside the observed la errorbar. This ensures 
that the effect of the different definition of ro is small. 
Finally, we note that equation M assumes that the dif- 
ferent data are uncorrelated. This is expected to be a 
fair assumption, especially for the X-ray data, which 
come from different catalogues. In the case of the opti- 
cal datasets, some larger correlation is probably present 
between the APM samples, due to the presence of com- 
mon objects in the two catalogues. However, we will 
show that we obtain almost identical conclusions about 
the confidence levels for the cosmological parameters 
considering the APM sample B only or the whole set of 
optical catalogues. This shows that the assumption of 
independence which underlies our likelihood analysis is 
not essential for the results. 

Finally, the best-fit cosmological parameters are ob- 
tained by maximizing £, i.e. by minimizing x^- The 95.4 
and 99.73 per cent confidence levels for the parameters 
are computed by finding the region corresponding to an 
increase A^2 with respect to the minimum value of x^- 
the exact value of A^2 depends on the given number of 
degrees of freedom z// . 

A similar analysis has been very recently performed 
by Robinson (2000), who used a simplified version of 
the above model and considered only APM data (all six 
subsamples). Wherever the comparison is possible, the 
results are qualitatively similar, but some quantitative 
differences are present. These differences can be ascribed 
to a different definition of the effective bias and, above 
all, to the absence of any account of the past-light cone 
effects in his analysis. 

We start our analysis by considering as a free pa- 
rameter f2oin only. We fix F = 0.2, which is in the 
range suggested by different works (see e.g. Peacock 
& Dodds 1996), and erg to reproduce the cluster abun- 
dance. For the normalization we adopt the fitting for- 
mula by Viana & Liddle (1999), who revised the Henry 
& Arnaud (1991) dataset of cluster X-ray temperatures 
and included a treatment of measurement errors. Their 
formula reads 



erg = 0.56r2on 



(7) 



where C = 0.34 and C = 0.47 for open and fiat mod- 
els, respectively. The claimed accuracy of this expres- 
sion is better than 3 per cent for r^om between 0.1 
and 1. The results of our maximum likelihood analy- 
sis obtained by using the complete set of catalogues 
are shown for flat models by heavy solid lines in Fig- 
ure hi The corresponding results for open models (not 
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99.73% 
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Figure 1. The variation of A 2 around the best-fitting value of the matter density parameter Oom for flat CDM models with 
shape parameter F = 0.2 and normalization reproducing the cluster abundance. The left panel refers to the results obtained 
using the optical dataset, while the right one refers to the dataset in the X-ray band. The heavy solid lines correspond to 
the results obtained using the complete cluster datascts, the light dotted ones refer to results obtained using only the largest 
catalogue (APM sample B and REFLEX, for optical and X-ray bands, respectively). Horizontal lines corresponding to the 95.4 
and 99.73 per cent confidence levels are also shown. 



shown here) are very similar. The left panel refer to 
the analysis of the optical catalogues. In this case the 
minimum of x^ is in the range fiom = 0.2 — 0.3 and 
we find r2omi$ 0.5 and r2omi$ 0.7 at 95.4 and 99.73 per 
cent levels, respectively. The results for the X-ray cat- 
alogues are shown in the right panel. Also in this case 
the most likely value for r2om is 0.2-0.3 but the allowed 
regions are much smaller, i.e. the X-ray data on the cor- 
relation length give tighter constraints: 0.2^ fiom^ 0.35 
and 0.2< ^om^ 0.45 at 95.4 and 99.73 per cent levels, 
respectively. This result is in qualitative agreement with 
our previous analyses, where we found that an Einstein- 
de Sitter model can be rejected because it predicts too 
low a correlation function for the RASSl and XBACs 
catalogues (Moscardini et al. 2000a, b). 

In order to understand which subsamples con- 
tribute to the constraints, we repeat the same analysis 
using in both the optical and X-ray cases the largest 
catalogues, i.e. APM sample B and REFLEX. The cor- 
responding curves are shown in the relative panels as 
light dotted lines. Unlike the case of the optical cata- 
logues where very small differences are found, the con- 
straints obtained using REFLEX only are less tight: 
0.2< fiom^; 0.5 and 0.2< Q.o^< 0.6 at 95.4 and 99.73 
per cent levels, respectively. This shows the importance 
of including in the analysis the subsamples at higher 
X-ray luminosities. 

In Figure H we show the results of our maximum 
likelihood analysis when we only fix the model normal- 
isation to reproduce the cluster abundance, leaving two 
free parameters: rJom and P. 

The upper panels refer to the analysis of the opti- 
cal catalogues. Confidence levels are shown separately 



for flat (left panels) and open (right panels) models. 
Here Uf = 2, i.e. A^2 = 6.17 and 11.8, for 95.4 and 
99.73 per cent levels, respectively. The allowed regions, 
at least for values of the density parameter f2om~ 0.4, 
appear to strongly depend only on P. The data sug- 
gest P ~ 0.15 ± 0.05 (errorbars are at 3(t level), both 
for open and flat models. At smaller r2om, the re- 
sulting values of P are larger, especially for models 
with non-zero cosmological constant. For instance, for 
f^Om = 0.2 we find P 



0.20 



+0.06 



0.24l°;J° for fiat models, and 
for open ones. The centre of the al- 



lowed region gives a relation between P and f2om , namely 
P = 0. 382-0. 797!r2om -fl.029!r2om -0.4780^ for the fiat 
models and P = Q.215-0.079r2om -0.0640^^+0. 069n{|,„ 
for the open models. 

The lower panels present the results obtained for 
the X-ray catalogues. Even if the allowed regions are 
qualitatively compatible with those displayed by the 
optical analysis, we notice some differences. First, the 
2- and 3cr regions are narrower, i.e. the X-ray data on 
the correlation length give tighter constraints on the 
cosmological parameters. Once again, the f2om depen- 
dence is small for high-density models: we find P in 
the range 0.1-0.15 (3a confidence levels) for f2om?^ 0.5, 
with no dependence on the presence of a cosmologi- 
cal constant. Again, models with low matter density 
favour larger values of P, between 0.2 and 0.3. The re- 
lations describing the centre of the allowed region are 
P = 0.487 - 1.342!r2om + 1.687nom - 0.73ng„, for the fiat 
models and P = 0. 394-0. 933r2om + 1.084ng„, -0.44f7g„, 
for the open models. 

Note that, combining the optical and X-ray data 
altogether, the allowed regions in the parameter space 
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Figure 2. Confidence contours (95.4 and 99.73 per cent confidence levels; light and heavy solid lines, respectively) of F and 
^Om for CDM models with normalization reproducing the cluster abundance. The left panels refer to flat cosmological models 
with varying cosmological constant Oqa = 1 ~ ^Om; the right ones to open models with vanishing SIqa- The upper row shows 
the results obtained using the optical dataset, while the lower one refers to the dataset in the X-ray band. The results obtained 
using only the largest subsamples (APM sample B for the optical case, REFLEX for the X-ray band one) are shown by light 
and heavy dotted lines (95.4 and 99.73 per cent confidence levels, respectively). 



(not reported in the Figure) are very similar to those 
obtained in the analysis of the X-ray catalogues only, 
showing once again that these data have larger discrim- 
inating power. 

Also in this case we checked the robustness of the 
resuhs when only the largest catalogues are used in the 
analysis. The results, presented by dotted lines, show 
that the confidence regions obtained using the APM 
sample B are very similar to those obtained using all 
three optical catalogues. On the contrary, when we use 
the REFLEX catalogue only, the resulting allowed re- 
gions, even if compatible with those obtained using also 



the catalogues at higher X-ray luminosity, appear to be 
larger. For instance, values of F as large as 0.4 cannot 
be excluded if the matter density parameter is small 

(fiom^ 0.2). 

In Figure H we present the constraints in the F — cts 
plane, once the values for fiom and fioA are fixed (i.e. 
Vf — 2). We consider an Einstein-de Sitter model, an 
open model with Slom = 0.3 and a fiat model, always 
with f2om = 0.3. Again the upper panels refer to the 
optical data. The confidence limits are quite similar for 
the three cosmological models, confirming once again 
the weak dependence on the density parameter. These 
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optical 




Figure 3. Confidence contours (95.4 and 99.73 per cent confidence levels; light and heavy solid lines, respectively) of F and 
erg for CDM models with (f2om, f^oA) = (li 0) (left panels), (0.3,0) (central panels) and (0.3,0.7) (right panels). The upper row 
shows the results obtained using the optical dataset, while the lower one refers to the X-ray dataset. The results obtained using 
only the largest subsamples (APM sample B for the optical case, REFLEX for the X-ray band one) are shown by light and 
heavy dotted lines (95.4 and 99.73 per cent confidence levels, respectively). 



results can be directly compared to the analysis made 
by Robinson (2000; its Figure 3); the agreement is good, 
even if at a given value of ag shghtly larger values of T 
are allowed in our analysis. 

The lower panels refer to the X-ray data and 
show the strongest constraints coming from these data. 
For example, for the Einstein-de Sitter model only a 
small region with values of F quite close to 0.1 and 
0.4^ (78^ 1.1 is allowed. The confidence limits obtained 
in the case of open and fiat models with rJom = 0.3 are 
similar and are consistent with the analysis of the op- 
tical data, but the allowed region is narrower: the 2g 
region has F in the range 0.14-0.22 and o^ between 0.6 



and 1.3 for the open model, while 0.16i; Fi; 0.22 and 
0.7^ (TgJ; 1.3 for the fiat model. 

Again the combination of optical and X-ray cat- 
alogues produces results almost indistinguishable from 
those obtained by the X-ray analysis only. The main 
difference is the further reduction of the allowed region 
in the case of the Einstein-de Sitter model: in this case 
we find F ~ 0.1 and 0.45< a8< 0.6. 

When the analysis is limited to the largest cata- 
logues, again no significant differences are found be- 
tween the confidence regions allowed by the APM sam- 
ple B and the complete set of optical catalogues. On the 
contrary, the results for the X-ray catalogues show that 
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the inclusion in the analysis of the shallower catalogues 
allows to reduce the confidence regions, excluding high 
values of F and as . 

In conclusions, the previous results show that the 
constraints coming from the X-ray datasets are in gen- 
eral tighter than those obtained from the optical data. 
This is also evident from the analysis limited to the 
largest corresponding catalogues (APM sample B and 
REFLEX). At first glance this is unexpected. In fact the 
errorbars of the correlation length ro for the APM sam- 
ple B are smaller than the REFLEX ones (see Table <^) 
and consequently their weights in the x^ estimates are 
larger. However, the galaxy clusters belonging to the 
X-ray datasets typically have a mass larger than the 
optical ones. In fact the density of the present optical 
clusters requires a minimum mass smaller than that ob- 
tained from the limiting fiuxes of the X-ray catalogues 
here considered. Since the bias factor entering in the 
correlation estimates is strongly dependent on the mass 
(see, for example, the discussion in Moscardini et al. 
2000b), this difference in mass produces a large spread 
of the ro values predicted for different cosmological mod- 
els, increasing the constraining ability of the results. For 
a similar reason, the inclusion in our analysis of other 
X-ray datasets besides REFLEX, even if with a smaller 
number of objects and with larger errorbars, helps in 
reducing the allowed regions in the parameter space. In 
fact these catalogues have a larger limiting flux than 
REFLEX and sample the cluster population at higher 
mass. 



5 PREDICTIONS FOR FUTURE SURVEYS 

5.1 The catalogues 

In this section we will use our model to predict the corre- 
lation length expected in possible future surveys both in 
the optical and X-ray bands. Thanks to their depth, in 
these surveys it would be possible to obtain information 
on the high-redshift behaviour of the cluster two-point 
function. For this reason we present our predictions di- 
viding the data in two different redshift bins {z < 0.3 
and z > 0.3), to allow a discussion of the redshift evo- 
lution of cluster clustering. 

In the optical band, a large improvement of our 
knowledge on the properties of the large-scale structure 
as traced by galaxy clusters will be obtained when the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) wiU 
be completed. This survey, which covers an area three 
times larger than the APM, will contain redshifts for ap- 
proximately one million galaxies. The expected number 
of galaxy clusters having at least 100 redshift measure- 
ments is approximately 1000. Moreover, the availability 
of 5 X 10^ galaxies in the photometric data (complete 
to a magnitude limit of r'=22) will allow the applica- 
tion of automated cluster-finding algorithms (e.g. the 
matched-filter approach), to extend the cluster cata- 
logue to higher redshifts. In its final form the SDSS 
is expected to be nearly as deep as two existing cat- 
alogues, the Palomar Distant Cluster Survey (PDCS; 
Postman et al. 1996; Holden et al. 1999) and the ESQ 



Imaging Survey (EIS; Olsen et al. 1999; Scodeggio et al. 
1999), which however are much smaller, covering only 
5.1 and 14.4 square degrees, respectively. These two cat- 
alogues have been built to find good distant candidates 
for successive more detailed observations and cannot be 
considered complete and well-defined samples suitable 
for statistical studies. Nevertheless, considering the ex- 
tended versions of these catalogues, the number density 
of candidates and the corresponding estimated redshift 
distributions are compatible. For these reasons we de- 
cided to use their characteristics to define the properties 
of a possible future cluster survey in the optical band. 
In particular, we use the number density measured for 
the EIS catalogue, which consists of 304 objects in the 
redshift range 0.2^ z'^ 1.3, with a median redshift of 
z ~ 0.5. The redshifts have been estimated by applying 
the matched-filter algorithm (Postman et al. 1996) and 
have an intrinsic uncertainty of at least Az = 0.1. A 
large effort began to validate these cluster candidates. 
Preliminary results (da Costa et al. 2000; Ramella et al. 
2000) show that more than 65 per cent of the studied 
candidates have strong evidence of being real physical 
associations. Moreover, first direct spectroscopic deter- 
minations of redshifts are in reasonable agreement with 
those derived from the matched-filter algorithm, with a 
possible systematic difference of Az ~ 0.1. We will com- 
pute the predictions of the clustering properties of this 
catalogue by fixing the minimum mass needed to repro- 
duce the EIS cluster density for z < 0.3 and z > 0.3. 
This choice automatically takes into account the pos- 
sibility that false candidates are included in the EIS 
catalogue by reducing the required minimum mass. As 
a consequence, the estimated correlation function will 
refer to the objects targeted as candidates and not to 
the real clusters which will be actually validated using 
spectroscopic techniques. The underlying assumption is 
that the failure of the technique used to select the candi- 
dates comes from the inclusion in the catalogue of high- 
luminosity systems with a mass just below the minimum 
mass of the true EIS clusters. 

In the X-ray band, the existing cluster clustering 
studies were expected to be largely overcome by the 
data that the ABRIXAS satellite (Triimper, Hasinger 
& Staubert 1998) was expected to collect, starting 
from mid 1999. Unluckily, problems with energy sup- 
ply caused the untimely loss of the satellite at the end 
of April 1999. In the plans, the ABRIXAS catalogue 
would have covered an area of 8.27 steradians up to a 
limiting fiux of Siim = 5 x 10~^^ erg cm~^ s~^ in the 
0.5-2 keV band. These characteristics have been used 
by Moscardini et al. (2000b) to make predictions of the 
cluster two-point correlation function. 

More recently, it has been proposed to use the very 
high sensitivity and good point-spread function of the 
XMM/Newton satellite, successfully launched in De- 
cember 1999, to build a very deep large-scale structure 
survey. The idea is to cover a region of 64 square degrees 
at high galactic latitude using 24 x 24 lOks XMM/EPIC 
pointings separated by 20 arcmin offsets. The expected 
limiting flux will be approximately Sum = 5 x 10~^^ erg 
cm~^ s~^ in the 0.5-2 keV band, which is 500 times 
more sensitive than the REFLEX one (see Pierre 2000 
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for more details). We predict the clustering properties of 
this sample using a constant sky coverage (in absence of 
real estimates) and assuming the previous limiting flux. 
Hereafter this sample will be called XMM/LSS. 

5.2 Results 

In Figure we show the predicted value of the correla- 
tion length for future surveys in the optical band (EIS) 
and in the X-ray band (XMM/LSS) as a function of 
the present matter density content flom- Here we con- 
sider only CDM models with as fixed to reproduce the 
cluster abundances (see equation 0) and we allow the F 
parameter to assume values from 0.1 to 0.5. The results 
are shown separately for clusters having z < 0.3 (upper 
panels) and z > 0.3 (lower panels). As a general result, 
we find that the presence of a cosmological constant in- 
creases the correlation length by a factor always smaller 
than 15 per cent. Once again we find that the depen- 
dence on the shape parameter F is strong: the higher 
F, the lower the predicted correlation length. On the 
contrary varying the matter density parameter, once F 
is fixed, changes ro only by a factor of at most 20 per 
cent. More precisely, we find a slight decrease of the cor- 
relation length with increasing Slom. All the catalogues, 
both in optical and in X-ray bands, display a positive 
redshift evolution of the clustering, i.e. the estimates of 
ro are larger for clusters at high redshifts. This result, 
which confirms a previous analysis by Moscardini et al. 
(2000b), is due to the large increase of the bias factor 
with redshift; this increase overcompensates the corre- 
sponding decrease of the dark matter correlation func- 
tion. In fact, considering high redshifts, galaxy clusters 
become rarer objects, connected to higher density fluc- 
tuations. 

The relatively low values of the predicted correla- 
tion length for the EIS catalogue show that its clus- 
ter number density corresponds to objects with a small 
mass, with some possible contamination coming from 
false candidates, which would reduce the amplitude of 
the expected clustering. 

Finally, we notice that the XMM/LSS sample has 
smaller ro than the ABRIXAS survey considered in our 
previous analysis (Moscardini et al. 2000b). Moreover 
the increase of clustering with redshift is less evident 
for the XMM/LSS sample than for the ABRIXAS cat- 
alogue. This is what we expect: when the limiting flux 
is decreased, a large number of small clusters enter in 
the catalogue, resulting in a smaller correlation function 
(see also Moscardini et al. 2000b). 

Figure n shows the predicted correlation length for 
the EIS catalogue when the geometry of the universe is 
fixed and the values of F and as are varied. We consider 
here the Einstein-de Sitter model and the open and flat 
models with f2om = 0.3. The small differences in the 
corresponding results confirm the slight dependence on 
the cosmology chosen. We find that, at a given erg, ro is 
a decreasing function of F, while, fixing F, ro is always 
an increasing function of as- 

A similar analysis has been performed using the 
characteristics of the XMM/LSS survey. Comparing the 
results, shown in Figure m, with the EIS ones, we can 



notice that the erg-dependence of ro is different: for very 
small normalisations, the correlation length is a decreas- 
ing function of as- This effect is due to the decrease of 
the effective bias, which is more rapid than the growth 
of the dark matter correlation function. 

It is interesting to note that, given a cosmological 
model, the predicted values of the correlation length for 
EIS and XMM/LSS are quite different, but the spread 
of ro obtained in different cosmological models is ex- 
pected to be of the same order of magnitude for the two 
cases. Unfortunately this does not allow to understand 
whether future X-ray or optical surveys will better con- 
strain cosmological parameters. 

Another important issue to be discussed is the con- 
straining ability of the clustering analysis with respect 
to the standard analysis based on cluster counts. 

In Figure W we show the expected number of galaxy 
clusters in the XMM/LSS survey as a function of the 
spectrum normalisation as- Results refer to the same 
cosmological models presented in the previous figures. 
We notice that these numbers are obtained directly from 
the model described above and do not pretend to re- 
produce the observed log A'^ — log S relation for clus- 
ters. To this aim it would be necessary to introduce an 
ad- hoc redshift evolution of the temperature-luminosity 
relation with one more parameter (see Section 3.2 and 
the discussion in Moscardini et al. 2000b). As already 
known, the results show a strong dependence on erg and 
a relatively weak dependence on F. For instance, using 
the values of the normalisation suggested by equation 
\A the change in the predicted number of clusters for 
0.1 < F < 0.3 (which includes the values suggested by 
a set of other observational data, see e.g. Peacock & 
Dodds 1996) is at most 50 per cent. If we analyse Fig- 
ure H we find that with the same assumptions (erg from 
cluster abundances and 0.1 < F < 0.3) the correlation 
length has a 100 per cent variation. This result allows us 
to conclude that clustering studies are a good comple- 
mentary tool to determine the cosmological parameters. 
In particular, once the normalisation erg is constrained 
by other observational data (i.e. cluster counts, cluster 
luminosity function or cosmic microwave background), 
it is quite useful in fixing the shape parameter F. 



6 CONCLUSIONS 

In this paper we discussed the model constraints that 
it is possible to infer from the analysis of the observa- 
tional data on the clustering of galaxy clusters. A set 
of 3 optical (coming from the APM and EDCC cata- 
logues) and 4 X-ray catalogues (RASSl Bright Sample, 
BCS, XBACs, REFLEX) has been considered. The the- 
oretical predictions for the different cosmological mod- 
els have been obtained by using a model which accounts 
for the clustering of observable objects in our past light- 
cone and for the redshift evolution of both the under- 
lying dark matter covariance function and the cluster 
bias factor. A linear treatment of redshift-space dis- 
tortions has been also included. In the case of X-ray 
selected clusters we followed the approach of Moscar- 
dini et al. (2000b), which makes use of theoretical and 
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Figure 4. The predicted value of the correlation length (in h^^ Mpc) for different future surveys in the optical band (EIS, left 
column) and in the X-ray band (XMM/LSS, left column) as a function of the present matter density content f2om. Results are 
presented for CDM models with normalization fitted to reproduce the cluster abundance. The different lines refer to different 
values of the F parameter: 0.1, 0.2, 0.3, 0.4 and 0.5, from top to bottom. Dotted lines: flat cosmological models (i.e. with 
non-zero cosmological constant); solid ones: open models with vanishing f2oA- The upper panels show the results obtained using 
clusters with z < 0.3, while the lower ones are for z > 0.3 clusters. 



empirical relations between mass, temperature and X- 
ray luminosity to convert the limiting flux of catalogues 
into a corresponding minimum mass for the dark mat- 
ter haloes hosting the clusters. In the optical band, we 
used the mean distance (i.e the number density) of the 
observed clusters to fix the minimum mass required by 
the model. 

Our theoretical predictions have been compared 
with the observed correlation lengths by means of a 
maximum-likelihood analysis. We considered cosmolog- 
ical models belonging to the cold dark matter class, de- 
fined by four parameters: the closure density in dark 



matter and vacuum energy (f2om and Qqa, respectively); 
the power-spectrum shape parameter F and normalisa- 
tion as- 

Our main results can be summarized as follows: 

• The constraints coming from X-ray and optical 
clusters are consistent but the former appear to be 
tighter. 

• If we fix the power-spectrum shape parameter F — 
0.2, in the range suggested by different observational 
datasets, and the power-spectrum normalisation to re- 
produce the cluster abundance, we obtain strong con- 
straints on the value of the matter density parameter, 
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Figure 5. The predicted values of the correlation length (in h^^ Mpc) for the EIS survey as a function of the value of the 
spectrum normalisation cf^. Results are presented for CDM models with (f2om, f^oA) = (li 0) (left panels), (0.3,0) (solid lines in 
the right panels) and (0.3,0.7) (dotted lines in the right panels). The different lines refer to different values of V: 0.1, 0.2, 0.3, 
0.4 and 0.5, from top to down. The upper panels shows the results obtained using clusters with z < 0.3, while the lower one 
are for z > 0.3 clusters. 



independently of the presence of a cosmological con- 
stant: nom < 0.5 and 0.2 < fiom < 0.35 at the 2a level, 
for the optical and X-ray data, respectively. 

• If we allow the shape parameter to vary, we find 
that the clustering properties of clusters are only weakly 
dependent on f^om and iloA. On the contrary, the results 
appear to be much more strongly sensitive to F. In fact, 
smaller shape parameters correspond to higher correla- 
tion lengths. 

• Considering models with the normalisation com- 
ing from the cluster abundance, we find that the cen- 
tre of the region allowed by the maximum-likelihood 
analysis of the optical data is described by the rela- 



0.797r2om + 1.029!:!L 



0A78Ql„ 



0.064r2L + 0.069r2g„ 



tion r = 0.382 

and r — 0.215 — 0.079r2om — u.uo'tssom 

for flat and open models, respectively. Considering the 

catalogues in the X-ray band, we find F = 0.487 — 

1.342nom + 1.687r2om - 0.73n(L for the flat models and 

F = 0.394 - 0.933!^om + 1.084nf,„, - O-Un^^ for the 

open models. 

• Using X-ray selected data only, wc find that for the 
Einstein-de Sitter model the value of F has to be quite 
close to 0.1 with 0.4^ ag'^ 1.1; for open and flat models 
with riom = 0.3 the 2cr region has 0.14< F< 0.22 and 
0.6< ag^ 1.3. 

We also used our model to make predictions on the 
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Figure 6. As Figure H but for the XMM/LSS survey. 



clustering properties of galaxy clusters expected in fu- 
ture surveys. We considered an optical catalogue with 
characteristics similar to the EIS project and a deep X- 
ray catalogue with the characteristics of the XMM/LSS 
survey. From this analysis we can conclude that: 

• Clusters at high redshifts are expected to have a 
larger correlation function than at low redshifts. 

• Again, predictions are almost insensitive to the 
presence of a cosmological constant while they are 
strongly dependent on the shape parameter. 

• The predicted clustering for the EIS catalogue is 
relatively small for all cosmological models suggesting 
that its cluster number density corresponds to objects 
with small mass, including some possible false candi- 
dates. 

• The correlation length for X-ray selected clusters is 



confirmed to depend on the limiting flux of the survey: 
the deeper the catalogue, the smaller ro. Moreover, the 
redshift evolution of clustering is less evident in deeper 
catalogues. 

In conclusion, our results show that the existing 
data on the clustering properties of clusters can be suc- 
cessfully used to put constraints on the cosmological 
parameters. The future availability of deeper surveys 
can increase the power of this approach, which can be 
considered complementary to the traditional study of 
cluster abundances. 
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Figure 7. The expected number of galaxy clusters with z < 0.3 (upper panels) and z > 0.3 (lower panels) in the XMM/LSS 
survey as a function of the value of the spectrum normalisation erg. Results are presented for CDM models with (f2om, f^oA) = 
(1,0) (left panels), (0.3,0) (solid lines in the right panels) and (0.3,0.7) (dotted lines in the right panels). The different lines 
refer to different values of the F parameter: 0.1, 0.2, 0.3, 0.4 and 0.5, from bottom to top. 
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